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"A method for determining the biological likelihood of candidate compositions or 

structures" 
Cross-Reference to Related Applications 

The present application claims priority from Provisional Patent Application No 
5 2003905362 filed on 1 October 2003, the contents of which is incorporated herein by 
referOTce, 

Field of the Invention 

This invention relates to a . method of detemiining the biological likelilaood of 
candidate compositions or structures, particularly glycans and Iheir derivatives also 
10 known as oligosaccharides. As used hereinafter the term glycan will include both 
glycans and glycan derivatives, unless otherwise indicated. 

Background of the Invention 

Glycans (sugar structures/oligosaccharides) are usually composed of varying 

15 numbers of less than a dozen biologically-occurring monosaccharides. When 
considered purely in terms of their masses there are usually only about 3-6 different 
mass-unique monosaccharides in a typical glycan structure. The most firequently 
encountered unique-mass monosaccharides are Hex (mass 162 Da; includes all hexose 
monosaccharides), HexNAc (mass 203 Da; includes all acetamidohexose 

20 monosaccharides), dHex (mass 146 Da; includes all deoxyhexose monosaccharides), 
Pent (mass 132 Da; includes all pentose monosaccharides), and NeuAc (mass 291 Da; 
N-acetylneuraminic (sialic) acid). There are several other biologically extant, though 
less-frequently encountered component monosaccharides, such as KDN, HexA, 
NeuGc. Other non-monosaccharide adducts such as sulfate (S; mass 79.97 Da), 

25 phosphate (P; mass 79.97 Da), methyl (14 Da), and acetyl (42 Da) are also occasionally 
observed on biologically-occurring oligosaccharides. 

It is often the case during the characterisation of biological molecules that a 
precise mass may be ascertained for each biological molecule but its composition and 
identity are unknown. Given a reasonably accurate mass, such as would normally be 

30 obtained by mass spectrometry, the monosaccharide composition of an unknown 
glycan can be theorised by determining, by computation, the set of monosaccharide 
compositions that are within a reasonable mass deviation (or tolerance) of the observed 
mass. This approach forms the basis of gly^comod 

(http://us.expasv.org/tools/glvcomod/) a publicly available research tool. The 

35 shortcomings of this tool, and of this purely theoretical approach by extension, is that a 
large number of compositions are returned for any mass of larger than moderate size. 
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and that the majority (90-99%) of these have little in common with known biologically 

extant compositions. 

The aim of the present invention is to attempt to alleviate some of the above 

described problems and to reduce the large number of irrelevant compositions returned 
5 by existing tools. 

Any discussion of documents, acts, materials, devices, articles or the like which 

has been included in the present specification is solely for the purpose of providing a 

context for ttie present invention. It is not to be taken as an adndssion that any or all of 

these matters form part of the prior art base or were common general knowledge in the 
1 0 field relevant to the present invention as it existed before the priority date of each claim 

of this application. 

Summary of the Invention 

In a first broad aspect, the preset invention incorporates statistical measures of 
1 5 biological relevance for the candidate compositions returned. 

Typically, biological relevance, expressed as a numerical score, or biological 
index, is determined by statistical comparison to an established reference set of known 
and fully characterised compositions, in the case of glycans a reference set such as the 
Glycosuite (http://www.glvcosuite.com) database. The biological index of any given 
20 composition may then be used as a basis for discarding biologically '^unlikely" 
compositions, as well as for ranking (sorting) of returned compositions by biological 
likeliness. 

Empirically, for glycans this allows between 90-99.9% of candidate 
compositions retumed by any given search to be discarded, whilst preserving and 
25 ranking the rraaaining, biologically likely compositions. 

In one aspect the present invention provides a method of determining the 
likelihood of a candidate composition comprising: 

selecting a reference group of known characterised compositions; 

establishing statistical characteristics relating to components of or other features 
30 of the known characterised composition; 

comparing the statistical characteristics of the known characterised 
compositions with corresponding components or features in the candidate compositions 
to establish a likelihood of those compositions occurring. 

35 More particularly, for glycans, in one aspect the present invention provides a 

method of characterising glycans comprising the steps of: 
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providing a search mass of a glycan whose composition is to be determined; 

generating a list of candidate glycans, typically including theoretical glycans, 
made up of components, including monosaccharides, whose total mass is within a 
predetermined tolerance of the search mass; 
5 selecting a reference group of known characterised glycan compositions of 

approximately similar mass to the search mass; 

establishing the mean and standard deviation of each component appearing in 
the reference group of the known characterised glycan compositions; 

for each candidate glycan composition calculating a partial score for each 
10 component in that glycan candidate based on the difference between the observed 
number of the component in the glycan candidate and the mean for that component in 
the reference group, divided by the standard deviation; 

combining the partial scores to provide an indication of the likelihood of that 
tiieoretical glycan candidate occurring. 
15 The candidate glycans will include structures which exist as well as theoretically 

possible glycan structures which are not known to exist. 

The predetermined tolerance of the search mass is within +/- 400Da, preferably 
+/- 200Da. 

The partial scores may be combined in any suitable maimer. One way, for 
20 example, is by multiplying the partial scores together. 

By the use of actual biological information, the present invention is able to 
discern biologically likely compositions from the vast majority of compositions of 
similar mass, but whose compositions are differ greatiy from known, biologically 
extant compositions. For example, for glycans where the publicly available web tool 
25 glycomod returns over 100 theoretical compositions for the mass 1300 Da +/- 0.5 Da, a 
tool embodjdng the present invention retums 2 biologically likely compositions and 
109 biologically unlikely compositions (which would normally be discarded). 

Although the main application of the present invention is to the delineation of 
biologically likely and imlikely sugar compositions for the purposes of sugar 
30 structure/composition elucidation, the generic methodology of using known biological 
data as a means to refine, interpret, and/or rank theoretical or empirical data may be 
used for many other applications. 
Brief Description of the Drawings 

Specific examples of the present invention will now be described by way of 
35 example only with reference to the accompanying drawings in which: 
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Figure 1 illustrates a computer running software embodying aspects of the 
present invention; 

Figure 2 shows the data entry page of a software based search tool embodying 
the invention; 

5 Figure 3 shows a results page of the search tool of Figure 2; 

Figure 4 shows the data entry page of the software based search tool for a mass 
that is known not to exist in the Glycosuite database; 

Figure 5 shows the results page for the search presented in Figure 4; 

10 Detailed Description of a Preferred Embodiment 

The present invention is implemented on a computer means running software 
carrying out the algorithms and process of the method. The computer is illustrated in 
Figure 1 and comprises a processor or CPU 100, a visual display screen 102, keyboard 
104, mouse 106 and printer 108. The computer is connected to a database 110 (known 

15 as GlycoSuite) by the internet a LAN or the like 

The first input to a search using a method to determine a glycan composition 
(the "search glycan") embodying the present invention, is a search mass (which is 
typically in Daltons). The search mass is typically an empirically determined mass of 
the "search glycaa" which is to be characterised determined by mass spectrometry or 

20 other means i.e. the mass of the search glycan whose composition is to be determined. 

A search mass tolerance (in Daltons) is also input Typically this will be a 
relatively small value depending on the expected accuracy of the empirically 
detemiined search mass and typically may be of tiie order of ± O.lDa. Also input is a 
"biological index" cut-oflF. The biological index is a measure of a theoretical glycan 

25 composition's likelihood and its derivation is explained in more detail below. The cut 
off is the value of that index above which candidate compositions are discarded as 
being too unlikely to occur in the real world. Also input is a ^^aximum composition" 
which indicates the maximimi allowable number of each monosaccharide in each 
theoretical glycan composition. By way of example, if it were known for a fact that the 

30 glycan to be characterised contained no sialic acid, the theoretical glycan compositions 
generated as potential matches for the search mass would also exclude sialic acid. This 
reduces the amount of computation required and improves speed and accuracy. In the 
system implementing the method, defaults would typically be provided for those inputs, 
except of course for the search mass. 



wo 2005/031343 



PCT/AU2004/001348 



5 

Other optional parameters may also be exposed to the user to further modify the 
performance of the search. The output of the composition search is a list of candidate 
glycan compositions, the majority of which will be theoretical compositions, i.e. 
possible structures but ones which are not known to be extant, whose mass is within, the 
5 search mass tolerance of the search mass, and whose biological index is less than the 
biological index cut-off. In theory one of those candidates matches the composition of 
the search glycan. The list of candidate glycan compositions may include naturally 
occurring glycan compositions. 

The composition search is performed as follows: 
10 Reference statistics for the given search mass are determined from the 

(Glycosuite) database. This process is described in more detail below. 

Monosaccharides are recursively recombined in varying numbers such that 
every possible combination of allowed monosaccharides is (treated. Compositions 
whose mass does not fall within the search mass tolerance are discarded, as are 
15 compositions for which the mmiber of any monosaccharide exceeds the maximum 
number of that monosaccharide specified by the "maximum composition". The result is 
a list of theoretical candidate glycan compositions. 

The biological index of candidate compositions is determined as described 
below. Compositions whose biological index does not satisfy the biological index cut- 
20 off are discarded. 

The remaining compositions are presented to the user in order of biological 
index. Typically the list will be short and may only include one or two candidates. 
This compares with the hundreds of candidates typically produced by Glycomod, each 
of which has to be individually reviewed and assessed. 
25 Calculation of Biological Index 

Inputs to the process are a composition, and a reference data set of known sugar 
compositions/structures. The reference set may be from any suitable database or data 
source such as Glycosuite. The output of the process is a numerical biological index. 

The determination of biological index for a given search glycan composition 
30 proceeds as follows: 

The mass of the composition is the search mass or may be determined by the 
sum of the residue masses of each monosaccharide/component in the composition. 

By reference to the reference set of known biological compositions, the mean 
and standard deviation of every monosaccharide/component in the database within an 
35 arbitrary mass range (eg: +/- 200 Da) of the mass of the composition is determined. 
Obtaining statistics from a range of masses around the given composition's mass is 
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necessary in order to obtain a sufficiently large sample size (preferably at least 100 
known compositions). In the case of the Glycosnite database of known sugar 
structures, a mass tolerance of 200 Da was empirically determined to be sufficient to 
provide in excess of 100 known compositions for search masses up to aroimd 3500. 

5 By way of example if the search mass were 1000 Da there may be 100 known 

glycans in the database whose mass is between 800 and 1200 Da. The mean and 
standard deviation of each of every monosaccharide/component appearing in those 
known glycans in the database is then determined. If we take HexNAc as an example 
we may iBnd that, on average, the 100 known glycans contain 3.3 HexNAc 

10 monosaccharides with a standard deviation of 2.3. This process is repeated to calculate 
the mean aiid standard deviation for each monosaccharide component Hex, dHex, pent 
et al, and each adduct in the known glycans, if adducts are being accounted for. 

For each candidate glycan composition "Partial scores" are then determined 
from the means and standard deviations calculated above. These are calculated for 

15 each monosaccharide in the given composition as the absolute value of the difference 
between the mean number of that monosaccharide in the reference set and the observed 
number of that monosaccharide in the theoretical candidate composition, divided by 
the standard deviation of that monosaccharide in compositions from the reference set. 
ie: 



20 



{mean — observed 



monosac 



stdev 

^'^^^ monosac 

where mearimonosac is the mean 'number of the given monosaccharide in the 
reference data set (Glycosuite); mearimonosac is the number of the given monosaccharide 
25 in the theoretical candidate composition; and stddevmonosac is fhe standard deviation of 
the given monosaccharide in the reference data set. 

By way of example if the theoretical glycan composition includes two HexNAc, 
three Hex and 1 NeuAc, the partial score for each of those three monosaccharides is 
calculated for that theoretical candidate glycan composition. Partial scores need not be 
30 calculated for monosaccharides which do not appear in the candidate theoretical glycan 
composition. 

In the event that the mearimonosac equals the mearimonosac for a particular glycan, 
flie system is arranged to give the partial score a minimum value of 0.01. 

Thus, the partial score of a monosaccharide is in fact the number of standard 
35 deviations the number of away from the mean that that monosaccharide is in the 
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candidate composition. In a normal distribution, approximately 68% of all data points 
lie within 1 standard deviation of the mean, '-93% within 2 standard deviations, over 
99% within 3. Assuming that the distributions of monosaccharide number for the mass 
range used to obtain the initial means and standard deviations for the given search mass 

5 are sufficiently close to normal, then partial scores of 3 or less for any monosaccharide 
indicate that the nxunber of those monosaccharides are within 99% of all compositions 
of similar mass in Glycosuite. 

Partial scores are then combined in some maimer to derive a single numeric 
score; this being the biological index. The actual mathematical derivation of the 

10 biological index may be arrived at using multiple means; different formulae exhibit 
subtlely different qualities in their sensitivity to large partial scores and other criteria. 
For this reason, biological index for the purposes of the present invention may be 
considered merely as a numerical value ttiat is representative o:^ and derived from, the 
magnitudes of the differences between a given composition and a population of known 

15 compositions of a similar mass. Presently, a biological index is calculated from partial 
scores as the product of all the partial scores from the candidate composition; ie: 
monosac^ 

monosacQ 

The Biological Index is adept at excluding very poor matches but at the same 
time if a candidate glycan composition has a very large (i.e. poor) partial score for one 

20 monosaccharide but low partial scores for the other monosaccharide components, the 
candidate may have an acceptably low Biological Index hence the system does not 
discard candidates which have only one poor partial score. 

The process of calculating the partial scores is carried out for each candidate 
glycan composition as discussed above. Compositions whose biological index does not 

25 satisfy the biological index cut-off are discarded. The remaining compositions are 
presented to the user in order of biological index. Typically the list will be short and 
may only include one or two remaining candidates. This compares with the hundreds 
of candidates typically produced by Glycomod, each of which has to be individually 
reviewed and assessed. 

30 The key element of the present invention is the use of biological data as a means 

to score the quality of theoretical data, in this case, sugar compositions. The actual 
manner in which a biological score/index is calculated is largely arbitrary; different 
formulae for calculating a biological index exhibit different characteristics with respect 
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to their tolerance to large compositional differences from the determined mean, and in 
their propensity to extrapolate the compositions present in the reference database. 

5 Example 1 

Figure 2 shows an entry page 200 of a software based search tool embodying 
aspects of the present invention. Three masses 202 are entered. In the example they 
are experimentally determined from bovine alpha-2-HS glycoprotein (fetuin). The 
entry page also allows the entry of tolerance 204, a biological index cut-off value 206 

1 0 and a maximum composition for various glycan components. 

Figure 3 shows the results page 300 for the search presented in Figure 2, listing 
all theoretical oligosaccharides compositions that are within the specified tolerance of 
each of the given search masses. Compositions that are present in the Glycosnite 
database are presented as hyperlinks to their relevant Glycosxiite record and are 

15 imderlined. Candidate compositions are presented in descending order of biological 
index. Compositions with a biological index of 2 or less are regarded as being very 
biologically likely. Those with a biological index of 2 to 10 are less likely but are 
considered to be sufficiently close to Glycosuite compositions of comparable mass. In 
the example, the search tool has identified the correct oligosaccharide composition and 

20 ranked it as the most biologically likely in each case. Biologically unlikely 
compositions (i.e. those having a biological index of greater than 10) are discarded and 
are not shown. 

Example 2 

25 Figure 4 shows the data entry page of the software based search tool for a mass 

400 that is known not to exist in the Glycosuite database. 

Figure 5 shows the results page 500 for the search presented in Figure 4. Three 
candidate oligosaccharide compositions are shown at 502, 504, 506, for the search mass 
2943 .4Da. The correct oligosaccharide composition is the third suggested composition 

30 506. All three suggested compositions have very low biological indexes indicating that 
all three compositions are deemed biologically likely, despite the fact that none of the 
three compositions is actually present in the Glycosuite database. This demonstrates 
that the searching tool can effectively extrapolate the Glycosuite database. Note that 
without using the biological index scoring and cut-off, there are 23 candidiate 

35 compositions that would satisfy the given (non-lenient) search criteria. With lenient 
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search parameters (i.e. unrestricted limits on monosaccharides), the number of 
theoretical compositions is over 100. 

5 

Example 3 

Reduced oligosaccharides released from a human respiratory mucin were 
detected by LC-MS. The compositions of the oligosaccharides were determined 
manually from MS/MS data, and the measured masses were submitted to searching by 

10 the software based search tool embodying the present invention( as known as 
"GlycoComp"). The masses were also subnMtted to GlycoMod ( a different glycan 
searching tool) for comparison. No limits were placed on the possible numbers of 
monosaccharides, except for methyl and acetyl groups, which were set to zero. Four of 
the compositions are preset in the GlycoSuite database, while the largest mass is not 

15 present. 
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20 Using a biological index cut-off value of 100 resulted in a reduction in the total 

number of likely compositions for each mass. This is particularly notable with large 
masses. The "GlycoComp" rank, sorted by biological index, also makes selection of the 
correct structure much easier, especially when compared with GlycoMod results, which 
are sorted only by mass. 

25 

Variations to the desaibed method are possible. For example, the partial score 
calculation as outlined above may be varied for improved accuracy. This improved 
monosaccharide partial score PartialScorem first calculates the number of standard 
deviations firom the mean for each monosaccharide: 

30 

StDevScorem= Abs(countnr meannO/stdeVm 
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This number of standard deviations from the mean is then converted to a 
probability, based on the normal distribution, for each monosaccharide: 



This improved partial score then represents the probability, based on similarity 
to the group of known oligosaccharides, that a particvilar number of monosaccharides 
are in fact present in a candidate oligosaccharide. The Biological Index is then the 
10 product of the partial scores for all monosaccharides. This method of calculating the 
Biological Index provides a combined total probability for the candidate composition. 
The biological index is then transformed to a larger number for easier interpretation. 



Partial scores are calculated based on the average characteristics of 
oligosaccharides present in GlycoSuiteDB within +/- 200 Da of the measured 
oligosaccharide mass. This molecular weight filter is based on the assumption that 
oligosaccharides with similar masses are more likely to have similar monosaccharide 
20 compositions. However, other filters based on biological source, disease state, or 
specific protein, could also be applied based on related assumptions — for example, that 
oligosaccharides from blood proteins, from homo sapiens, with cancer may have 
similar monosaccharide compositions. These filters could be user-defined in various 
combinations. 

25 Filters could also incliide more processed measures, such as the ratios of 

monosaccharides in an oligosaccharide (eg fiicose:hexose, HexNAciall 
monosaccharides). Rather than averaging values after a particular filter, equations of 
best fit could also be used to calculate partial scores (eg the ratio of fiicose:hexose as a 
function of molecular weight). These could also be used in various combinations (eg 

30 the ratio of fucose:hexose is a function of molecular weight, with a filter for only 
oligosaccharides from mammals). Variance associated with these lines of best fit (eg R 
values) could be used to calculate partial scores. 

Since oligosaccharides from the same biological source are likely to be 
structurally and compositionally related, a nested biological index (NBI) could be used 



5 PartialScorem= 



— iStDevScore„)^ 

e ^ 



-Jin X stdev, 



BiologicalInd< 
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to determine how similar the best ranked compositions, based on biological index, are 
to each other. This might be particvilarly useful to suggest cases in which a composition 
has a high biological index, due to similarity to oligosaccharides in the GlycoSuite 
database, but is quite different to the compositions retumed for other search masses 
5 from the same sample. Since experimental oligosaccharides from one sample are likely 
to have a wide range of masses, a mass filter is unlikely to be useful in this case. 
However, other measures, such as ratios of monosaccharides, may be less dependent on 
molecular weight. The best ranked compositions for each mass could constitute a 
'database', to which each individual composition could be compared. This comparison 

10 could be made by a NBI, or other analysis, such as principal component analysis 
(PCA). If there were any matched compositions (by biological index) which were 
significantly different to the other matches, this could indicate that they were not the 
correct composition, based on the biological characteristics of the sample. 
Compositions ranked lower by biological index could then be included in place of the 

15 putatively incorrect matches, and NBIs again calculated. This could proceed until there 
were minimal matches that were different to other retumed matches, (eg until the 
standard deviations of each monosaccharide for all matched compositions were 
minimised). 

A further extension to the calculation of biological index would be to 
20 incorporate estabUshed rules into the calculation of probabilities. Calculating the 

percentage of structures within a window that conform to a particular rale and 

converting this to a percentage could achieve tiiis. 

To automatically account for glycopeptide mass additions to an oligosaccharide, 

the software based search tool could either digest a known glycopeptide with a 
25 protease, or add the amino acid residue masses to the software based search tool's 

alphabet in order to predict possible amino acid compositions in addition to 

oligosaccharide compositions. The number of possible amino acid compositions may 

be limited by requiring the presence of certain amino acid motifs that sugars can attach 

to. 

30 Although the present invention as described above is concerned with the use of 

known sugar structures/compositions as a means to discern/elucidate monosaccharide 
composition given only a mass, the concept could be extended to other compositions 
and to the use of other stractural characteristics, for example linkage and branching, as 
reference data for determining and/or ascertaining the quality of complete sugar 

35 structures for other investigative techniques, such as glycan fragment mass 
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fingerprinting (see the applicant's co-pending provisional patent application No 
2003902907, the entire contents of which are incorporated herein by reference). 

It will be appreciated by persons skilled in the art that numerous variations 
and/or modifications may be made to the invention as shown in the specific 
5 embodiments without departing firom the spirit or scope of the invention as broadly 
described. The present embodiments are, therefore, to be considered in all respects as 
illustrative and not restrictive. 



